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THREE CONSECUTIVE ALMOST SQUARES 


JEREMY ROUSE AND YILIN YANG 


Abstract. Given a positive integer n, we let sfp(n) denote the squarefree part of n. We 
determine all positive integers n for which max{sfp(n), sfp(n + 1), sfp(n + 2)} < 150 by relating 
the problem to finding integral points on elliptic curves. We also prove that there are infinitely 
many n for which 

max{sfp(n),sfp(n + l),sfp(n + 2)} < 


1. Introduction 

The positive integers 48, 49 and 50 are consecutive, and are “almost squares”, namely 48 = 3-4^, 
49 = 7^ and 50 = 2 • 5^. Does this phenomenon ever occur again? That is, is there a positive 
integer n > 48 for which 

n = 

n + 1 =y‘^ 

n + 2 = 2z^ 

has an integer solution? The answer is no. One perspective on this is given by Cohen in Section 
12.8.2 of [5], where this problem is solved using linear forms in logarithms. Another approach 
is to recognize that the pair of equations 2z^ — = 1 dehne an intersection of two 

quadrics in and hence dehne a curve of genus 1. Siegel proved that there are only hnitely 
many integer points on a genus 1 curve, and hence the system of equations above has only 
hnitely many solutions. This method of solving simultaneous Pell equations is considered in 

BI]. 

We consider more generally the problem of, given integers a, b and c, hnding integers n for which 
n = ax"^, n + 1 = by‘^ and n + 2 = cz^. Multiplying these equations gives n{n + l)(n + 2) = 
abc{xyzY and is related to the problem of hnding consecutive integers that multiply to an 
“almost” perfect square. There is an extensive literature on related problems. On one hand, 
Erdos and Selfridge proved in [7] that no product of consecutive integers is a perfect power (a 
result generalized to arithmetic progressions by others, see |1] and |I2]). In another direction, 
Gydry proved HH that the equation n{n + 1) ■ ■ • (n + fc — 1) = bx^ has no integer solutions with 
n > 0 if the greatest prime factor of b does not exceed k. 
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For a positive integer n, define sfp(?7,) to be the “squarefree part” of n, the smallest positive 
integer a so that a|n and - is a perfect square. In this paper, we consider positive integers n 
for which sfp(n), sfp(n + 1) and sfp(n + 2) are all small. Our first result is a classification of 
values of n for which sfp(n), sfp(n + 1) and sfp(n + 2) are all < 150. We say that such an n is 
“non-trivial” if sfp(n) < n, sfp(n + 1) < n + 1 and sfp(n + 2) < n + 2. 

Theorem 1. There are exactly 25 non-trivial n for which sfp(n) < 150, sfp(n + 1) < 150 and 
sfp(n + 2) < 150, the largest of which is n = 9841094. 

Remark. A table of all 25 values of n is given in Section\^ 

After this manuscript was complete, the authors were informed that the essence of this result 
follows from work of Michael Bennett and Gary Walsh (see H). Given positive integers a, b 
and c, a positive integer solution to n = n + 1 = by"^, n+2 = cz"^ yields (n + 1)^ —n(n + 2) = 
b'^y^ — {ac){xzY = 1- Bennett and Walsh prove that the system — dX"^ = 1 has at most 
one solution in positive integers. Moreover, if T + U\fd is a fundamental solution of Pell’s 
equation — dY'^ = 1 define for fc > 1, + Uk\fd = (T + Uy/d)’^. Then = bx'^ for some 

integer x G Z for at most one k, and if any such k exists, then k is the smallest positive integer 
for which is a multiple of b. Bennett and Walsh’s approach uses linear forms in logarithms, 
while our approach relies heavily on the theory of elliptic curves. In particular, if C is the curve 
defined by the two equations by^ — ax^ = 1, cz^ — by‘^ = 1, then the Jacobian of C is isomorphic 
to E : y"^ = x^ — {abcYx. 

We rule out many of the 778688 candidates for the triple (a, b, c) by checking to see whether 
there are integer solutions to each of the three equations by"^ — ax^ = 1, cz^ — 6^ = 1 and 
cz"^ — ax^ = 2. We also test C for local solvability, and use Tunnell’s theorem (see ra) 
to determine if the rank of E is positive. Finally, we use the surprising property that the 
natural map from C to E sends an integral solution on C to an integral point on E. It suffices 
therefore to compute all the integral points on E (using Magma |2]), which requires computing 
generators of the Mordell-Weil group. In many cases this is straightforward, but a number of 
cases require more involved methods (12-descent, computing the analytic rank, and the use of 
Heegner points). 

Given the existence of large solutions, it is natural to ask how large 

max{sfp(?7,), sfp(?7, -|- l),sfp(n -|- 2)} can be as a function of n. This is the subject of our next 
result. 

Theorem 2. There are infinitely many positive integers n for which 

max{sfp(n), sfp(n -|- 1), sfp(u -|- 2)} < 

Remark. The following heuristic suggests that the exponent 1/3 above is optimal. A partial 
summation argument shows that the number of positive integers n < x with sfp(n) < z is 
-|- 0{^/x\og{z)). Assuming that n, n 1 and n 2 are “random” integers it follows 
that sfp(n), sfp(n -|- 1) and sfp(n -|- 2) are all < n" with probability about Therefore, 

the “expected” number of n for which max{sfp(?7,), sfp(n -|- l),sfp(?7, -|- 2)} < u" is infinite if 
3(a — l)/2>—1 and finite otherwise. 
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2. Background 


We denote by Qp the field of p-adic numbers. A necessary condition for a curve A/Q to have 
a rational solution is for it to have such a solution in Qp for all primes p. 

If d is an integer which is not a perfect square, let Z[\/d] = {a + b\fd : a,b ^ Z}. This is a (not 
necessarily maximal) order in the quadratic field Q[\/d]. Let N : Z[\/d] —)■ Z be the norm map 
given by A(a + bVd) = of — db'^. 

We now describe some background about elliptic curves. For our purposes an elliptic curve is 
a curve of the shape 

E : y^ = ax b 

where a,b E Q. Let E{Q) be the set of pairs {x, y) of rationals numbers that solve the equation, 
together with the “point at infinity”. This set has a binary operation on it: given P, Q in P(Q), 
the line L through P and Q intersects P in a third point R = (t, y). The point P + Q is defined 
to be {x, —y). This binary operation endows P(Q) with the structure of an abelian group. 

The Mordell-Weil theorem (see mi. Theorem VIIL4.1, for example) states the following. 
Theorem 3. For any elliptic curve P/Q, the group P(Q) is finitely generated. More, precisely, 

P(Q) ^ P(Q)tors X 

where E{Q)tors is the (finite) torsion subgroup. 


There is in general no algorithm which is proven to compute the rank of P, (see Section 3 
of Rubin and Silverberg’s paper |T6]) but there are a number of procedures which work well 
in practice for relatively simple curves P. We will reduce the problem of solving n = ax^, 
n + 1 = by'^ and n + 2 = cz^ to finding points {X, Y) on the curve 

E:Y^ =X^- (abcfX, 


with X, y G Z. A theorem of Siegel (see Theorem IX.4.3 of mi) states that there are only 
finitely many points in P(Q) with both coordinates integral. There are effective and practical 
algorithms (see [19], [9], [T8| and [22]) to determine the set of integral points, provided the rank 
r can be computed and a set of generators for P(Q) found. Given a point P = {x,y) G P(Q), 
the “naive height” of P is defined by writing x = | with a, 6 G Z with gcd(a, b) = 1 and defining 
h{P) = logmax{|a|, |6|}. The “canonical height” of P is defined to be 


h{P) 


lim 

n—^oo 


hifMP) 

4 ?! 
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The L-function of T^/Q is defined to be 

an{E) 


HE,s) = J^ 


JJ (1 - a^(E)p ‘ + £{p)r/ ^‘) ‘ 


n=l 


n° 


where ap{E) = p + 1 — jj^E{¥p), and 

1 if p does not divide the conductor of E^ 


e{p) = 


0 otherwise. 


The Birch and Swinnerton-Dyer conjecture states that ords=iL(ii^, s) = rank(ii^(Q)), and more¬ 
over that 

L{E,s) f^i?(i?/Q)m(E/Q)n cp 

hm -— =-—- - —. 

(s — !)'■ |Tp 

In our case, hi is twice the real period, R{E/Q) is the regulator of E{Q) computed using 
the function h above, the Cp are the Tamagawa numbers, and III(T^/Q) is the Shafarevich- 
Tate group. The completed L-function A(E,s) = N^^'^{2 ti)~^T{s)L(E, s) satishes the function 
equation A{E,s) = weA{E,2 — s), where we is the root number of E. Note that we = 1 
implies that ords=iL(L^, s) is even, and we = —1 if ords=iL(L^, s) is odd. 

The best partial result in the direction of the Birch and Swinnerton-Dyer conjecture is the 
following. 

Theorem 4 (Gross-Zagier [10], Kolyvagin [14], et ah). Suppose that E/Q is an elliptic curve 
and oYds=iL{E, s) = 0 or 1. Then, oYds=iL{E,s) = rank(i?(Q)). 

The work of Bump-Friedberg-Hoffstein [1] or Murty-Murty [H] is necessary to remove a con¬ 
dition imposed in the work of Gross-Zagier and Kolyvagin. 


3. Proof of Theorem [2] 

Proof. Motivated by the observation that 8388223 = 127 • 257^ and 8388225 = 129 • 255^, we 
find a parametric family of solutions where n = a ■ and n -|- 2 = (a -|- 2) ■ (6 — 2)^. 

If we write (4 -|- \/T3)(649 -|- 180\/l3)”^ = Xm + 2/m\/l3 where Xm and ym are integers, then 
— I3y^ = 3 for all m > 0. We have that Xq = 4, xi = 4936, y^ = 1, yi = 1369 and 

Xm 1298Xm—1 Xm—2t ym 1298?/^—1 |/m—2- 

It is easy to see that Xm is periodic modulo 32 with period 8 and from this it follows that 
Xgm+y = 0 (mod 32) for all m > 0. Set Om = x^m+i and n = 4a^ — 3am — 1- Then we have 

n iflm l)(2nm T 1) > 

n + l = am(4a^ - 3) = l3amylm+7, and 

n -\- 2 = (om + l)(2am — 1)^- 
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Since 16|am, max{sfp(n), sfp(?7,+l), sfp(?7,+2)} < max{am —1, Om+l} = Om+l < □ 

Remark. The polynomial — 3x used in the proof above is the Chebyshev polynomial T^{x). 
This explains why T^i^x) — 1 and T^i^x) + 1 both have a double zero. 

4. Proof of Theorem [U 

Since there are 92 squarefree integers < 150, there are 778688 = 92^ possibilities for the triple 
(a, b, c). We hrst test four things before searching for integral points on B : y^ = x^ — (abc)^x. 
Suppose that n = n + 1 = by'^ and n + 2 = cz^ is an integral solution to the system of 
equations 

( 1 ) by^ — ax'^ = 1 , 

( 2 ) -by'^ = 1, 

(3) cz"^ — ax^ = 2. 

4.1. Greatest common divisor conditions. If {x,y,z) is an integer solution to ([1]), ([2]) and 
([3]), then gcd(a, b) = 1, gcd(6, c) = 1 and gcd(a, c) = 1 or 2. This reduces the number of triples 
to consider to 425639. 

4.2. Norm equations. We next check that each of the equations by'^ — ax"^ = 1, cz"^ — by'^ = 1, 
and cz"^ — ax"^ = 2 has an integer solution. The equation by'^ — ax^ = 1 has an integer solution 
if and only if — abx^ = b has an integer solution. This equation has a solution if and only if 
Z[\/^] has an element of norm b. We use Magma’s routine NormEquation to test this. After 
these tests have been made, there are 2188 possibilities for (a, b, c) that remain. 

4.3. Local solvability of C. Let G C be the curve dehned by the two equations by'^ — ax^ = 
cz"^ — by^ = w"^. This curve must have a rational point on it in order for there to be a 

non-trivial solution. We check whether C has points {x \ y z \ w) in for all primes p 
dividing 2abc. This is done via Magma’s IsLocallySolvable routine. This eliminates 244 
possibilities. 

4.4. Rank of the elliptic curve E. Let E \ y'^ = x^ — {abcYx. Dehne a map M from non¬ 
trivial solutions to G, represented in the form n = ax"^, n-|- 1 = by^ and n + 2 = cz^ to E, given 
by 

M{x, y, z, n) = ((n -|- l)(a6c), (abcYxyz). 

If n > 0, then xyz > 0 and so M{x,y, z,n) G E{Q) is an integral point, and one with a 
non-zero y-coordinate. We now use the following results about the family of elliptic curves E 
we consider. 

Definition 5. A natural number N is called congruent if there exists a right triangle with all 
three sides rational and area N. 
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Consider the elliptic curve E over Q given by: 

E{Q) : 7/2 = ^3 _ 

We have the following result. 

Theorem 6 (Proposition 1.9.18 of M). The number N is congruent if and only if the rank of 
E is positive. 


Proposition 1.9.17 of [13] implies that the only points of hnite order on E are (0, 0), (±A^, 0) and 
the point at inhnity. If for some n > 0 we have n = n + 1 = hy"^ and n + 2 = cz"^, then the 
point M{x,y,z,n) on E : — {abc)‘^X is, according to the above, a non-torsion point, 

and hence E has positive rank, and consequently, by Theorem 16] abc is a congruent number. 


Theorem 7 (Tunnell 1983). If N is squarefree and odd, and 

# {x, y,z E 'L\N = 2x^ + |/^ + 32z‘^] ^ -# {x, y,z E 'L\N = 2x^ + + 8;^^} , 

then N is not congruent. If N is squarefree and even, and 
# |x,|/, 2 : e Z : 
then N is not congruent. 


N 0 9 9 

— = Ax^ + y^ + 32^2 


1 


N 


7 ^ -# < x,y, z E ■. — = 4x + y + 8z- 


To use Tunnell’s theorem, we compute the generating function for the number of representations 
of N by 2x2 + y‘^ + 8z^ as 



2x^-l-y^-h8z^ _ 


x,y,zGZ 






as well as the representations of N by 2x^ + y^ + 32z^, 4x2 + y^ + 8z^ and 4x2 + y^ + 32z^. 
Actually, for our purposes, it suffices to compute 

X Y Z 




^2x^+y^+8z^ 


x=-Xy=-Y z=-Z 

where X = [^150V2j, Y = and Z = [^ISOVSJ. 


The use of Tunnell’s theorem rules out 530 of the 1944 remaining cases. In the case that the 
hypothesis of Tunnell’s theorem is false, the Birch and Swinnerton-Dyer conjecture predicts 
that £'(Q) does have positive rank. In this case, we proceed to the next step. 


4.5. Computing integral points. Once the program determines that the elliptic curve i7(Q) : 
y2 = — {abcYX has positive rank, we let Magma attempt to compute the integral points on 

the curve, using the routine IntegralPoints, which is based on the method developed in [9] 
and [T9|. If this routine does not succeed within 15 minutes, which is sufficient time to perform 
a 4-descent, we abort the computation. 
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In 1377 of the 1414 cases that remain, Magma is able to determine the Mordell-Weil group and 
determine all of the integral points within 15 minutes. Once the integral points are determined, 
we check to see if they are in the image of the map M{x,y,z,n) and if so, whether they 
correspond to a non-trivial solution. 

There are 37 more difficult cases that remain, and in each case we are able to use other methods 
to determine the Mordell-Weil group. A table of these cases and the generators of the Mordell- 
Weil group are given at the page http://users.wfu.edu/rouseja/MWgens.html, 

Of the 37 cases, there are four cases with root number —1 and rank < 3 for which one point 
of infinite order has low height. In these, we numerically compute L'{E, 1) and show that it 
is nonzero. Theorem 0] proves the rank is one in these cases. One of these cases is a = 139, 
b = 89 and c = 109. In this case, E has conductor 5.82 • 10^^ and computing L'[E, 1) takes 
about 5 hours. 

In the remaining 33 cases, rather than use the 8-descent (the default if the IntegralPoints 
command were to keep running), we use the 12-descent routines in Magma (due to Tom Fisher 
[H]) to hnd points of large height. The curve with a = 137, b = 109 and c = 101 has a 
generator with canonical height 1234, which is found after searching on the 12-covers for about 
90 minutes. This approach is successful in all but one case. 

4.6. Heegner points. All but one of the 778688 original cases are handled by the methods 
of the previous sections. The remaining case is a = 67, b = 131 and c = 109. The curve E 
has root number —1, rank < 1 and conductor 2.9 ■ 10^^. A long computation shows that 
L'{E, 1) 72.604. This suggests, assuming the Birch and Swinnerton-Dyer conjecture is true 

and in(i7/Q) is trivial, that a generator of the Mordell-Weil group has canonical height about 
1692. Searching for points on the 12-covers up to a height of 3^^ • 10® ~ 10^® does not succeed 
in finding points. 

For this hnal case, we use the method of Elkies described in [6]. This is a variant of the usual 
Heegner point method and is quite fast on quadratic twists of curves with low conductor. The 
modular curve Xo(32) parametrizes pairs {E,C), where E is an elliptic curve, and C C E is 
a cyclic subgroup of order 32. It is well-known that Xo(32) is isomorphic to y‘^ = + 4a;, 

which in turn has a degree 2 map to y"^ = x^ — x. Finding a rational point on y'^ = x^ — D^x 
is equivalent to Ending a point on y'^ = x^ — x with x and y/\/—D both rational. If E is an 
elliptic curve with complex multiplication whose endomorphism ring O contains an ideal I so 
that Ojl = Z/32Z, this naturally gives rise to a point on Xo(32), which is defined in the ring 
class field of O. Taking the trace of such a point (in the Mordell-Weil group) gives rise to a 
point on Xo(32) over an imaginary quadratic field. 

In [6], Elkies gives a procedure for constructing rational points on y‘^ = x^ — D‘^x with D = 7 
(mod 8). We must handle one case with D = 5 (mod 8). We let O = Z[4\/— //]. This ring 
has an ideal / = (32,4 -|- 4\/— 77) and Ojl = Z/32Z. We take a representative ideal J for 
each element of the class group of O. Thinking of J as a lattice, C/IJ is an elliptic curve, and 
J//J is a subgroup of order 32 on the curve. Thus, {C/IJ, I/U) is a Heegner point on Xo(32). 
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Via the isomorphism with E ■. + 4a;, we obtain a collection of points. Adding these 

points together on E gives a point dehned in Q{y/—D), and eventnally a rational point on 
— D'^x. Rnnning this compntation in several simpler cases snggests that the resnlting 
point P on y^ = — D^x satisfies mQ = P for some rational point Q, where m is eqnal to 

the nnmber of divisors of D. Any point Q must have an x-coordinate of the form x = 
where e is a divisor of 2D. We use this method to minimize the amount of decimal precision 
needed to compute P. 

For D = 67-131-109, the class number of O is 3712. Given that we expect a point with canonical 
height 1692, we compute the Heegner points and their trace in the manner described in [6] using 
850 digits of decimal precision. The result is a point Q with canonical height 1692.698, and the 
computation takes about 18 minutes. This point is a generator of the Mordell-Weil group of 
E : y'^ = x^ — D‘^x. Once the Mordell-Weil group is found, we check that there are no integral 
points on E. This completes the proof of Theorem [T] 

5. Table of n 

The following is a table of all 25 positive integers n with sfp(fc) < max{/c, 150} for k = n, n + 1 
and n -|- 2. 


n 

sfp(n) 

sfp(n + 1) 

sfp(n + 

48 

3 

1 

2 

98 

2 

11 

1 

124 

31 

5 

14 

242 

2 

3 

61 

243 

3 

61 

5 

342 

38 

7 

86 

350 

14 

39 

22 

423 

47 

106 

17 

475 

19 

119 

53 

548 

137 

61 

22 

845 

5 

94 

7 

846 

94 

7 

53 

1024 

1 

41 

114 

1375 

55 

86 

17 

1519 

31 

95 

1 

1680 

105 

1 

2 

3724 

19 

149 

46 

9800 

2 

1 

58 

31211 

59 

3 

13 

32798 

62 

39 

82 

118579 

19 

5 

141 

629693 

53 

46 

55 

1294298 

122 

19 

7 

8388223 

127 

26 

129 

9841094 

134 

55 

34 
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